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A TECHNIQUE FOR DESIGNING ACTIVE CONTROL SYSTEMS 
FOR ASTRONOMICAL-TELESCOPE MIRRORS 

By W. E. Howell and J. F. Creedon 
Langley Research Center 

SUMMARY 

* 

This paper considers the problem of designing a control system to achieve and 
maintain the required surface accuracy of the primary mirror of a large space telescope. 

Control over the mirror surface is obtained through the application of a corrective force 
distribution by actuators located on the rear surface of the mirror. The design proce- 
dure is an extension of a modal control technique developed for distributed parameter 
plants with known eigenfunctions to include plants whose eigenfunctions must be approx- 
imated by numeric-! techniques. Instructions are given for constructing the mathemat- 
ical model of the system, and a design procedure is developed for use with typical numer- 
ical data in selecting the number and location of the actuators. 

Two techniques for treating disturbances to the plant are discussed. These two 
techniques, which treat the errors as deterministic and uncorrelated, respectively, are 
examined from the standpoints of sensitivity to various mirror errors, determining the 
number of actuators required, and means of finding the best locations. For the deter- 
ministic case it was found that the "best" actuator locations (those locations which will 
minimize the steady-state error) aie very sensitive to the error distribution. In addition, 
these locations can presently be found only by exhaustive searches of all possible actuator 
locations, and the number of actuators required for a specific mirror and specific error 
can only be estimated after much computer time is used. In practice the error distribu- i 

tion over the mirror surface would be expected to change with the telescope attitude rela- 
tive to the sun. Also, the exact nature of the mirror errors will N» ti.ne varying and will 
not, in any case, be known very precisely. For these reasons it is not recommended that 
the errors be treated deterministically. In addition, when the errors at any particular 
time are treated as uncorrelated random variables, the actuator locations are much less 
sensitive to specific variations in error distribution, an estimate of the number of actua- 
tors required to produce a desired reduction in figure error can easily be made, and loca- 
tions which will yield results near the estimated figure accuracy can be found in a reason- 
able manner. Thus, at present this technique is preferred even though it requires more 
actuators than the deterministic method for a specific assumed error. 




Several numerical examples are presented for a 76-cm-diameter (30-inch), thin 
spherical mirror and the computer program to implement the design procedure is given 
in an appendix. The results include a comparison of the modal control law and an opti- 
mal (least-squares) control law. The results of this comparison indicate that not much 
performance is to be gained by the added complexity of this optimum control law. 

INTRODUCTION 


One of the mcst fundamental problems associated with orbiting a large, diffraction- 
limited telescope of the size and type discussed in reference 1 i3 that of manufacturing, 
figuring, and maintaining the figure of the large primary mirror. Many factors such as 
initial figuring errors, the change from lg to Og, and changing temperature gradients on 
the mirror while in orbit make conventional techniques of figuring and supporting tele- 
scope mirrors unsuitable. An alternate approach has been developed in which the mirror 
is actively controlled by first sensing figure errors on the primary mirror and then null- 
ing them by properly deforming the mirror. This technique (ref. 2) has been successfully 
applied to a 76-cm-diameter (30-inch), 1.27-cm-thick ^-inchj mirror with an initial 

error of 1/2 wavelength rms (X = 0.6328 pm). By using 56 actuators, the mirror was 
controlled to a figure accuracy of better than /Id. 

Since the time of the investigation reported in reference 2, consideration has been 
given to the application of a modal control technique to a class of mirrors. The modal 
control technique represents the plant to be controlled in terms of the eigenvalues and 
eigenfunctions of the linear differential operator which describes the behavior of the 
plant. In reference 3, this technique was applied to distributed parameter plants whose 
eigenvalues and eigenfunctions could be obtained in closed form. In many practical 
examples, however, the required eigenfunctions are not available. For the mirror, for 
example, the restrictions of practical mounts (boundary conditions) and the existence of 
holes in the center of the primary mirror used in Cassegrain telescopes preclude obtain- 
ing the required eigenfunctions in closed form. Estimates of these functions must be 
obtained via numerical approximation techniques. The purpose of this paper is to set 
forth for such a system a design procedure based on the use of the modal control law 
described in reference 3. First, the modal control concept is explained, the control sys- 
tem described, and the analysis procedure set forth. Certain specific details, such as 
accounting for the pad effects and the treatment of initial or expected error, are then 
covered. Numerical data for the 76-cm-diameter (30-inch) mirror are given and exam- 
ples are presented. Appendix A contains a listing of the computer program; appendixes B, 
C, and D, eigenvector listings and diagrams and eigenvalue listings for the mirror; appen- 
dixes E and F, several examples of actuator placement and resultant mirror errors. 
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SYMBOLS 


Values in the body of the paper are given both in SI Units and U.S. Customary Units. 

The measurements and calculations were made in the U.S. Customary Units. The values 

in the appendixes are in U.S. Customary Units and are consistent with the program in 

appendix A. 

A area of mirror 

aj ith modal coefticient which expands P in terms of the mode shapes (see 

eq. (35)) 

a N N x 1 vector of coefficients of the force distribution in the modal domain 

which corresponds to the controlled modes 

a** R x 1 vector of coefficients of the force distribution on the mirror; these 

coefficients arise from the action of the N actuators 

C M x 1 vector which is the sum of the control system displacements and the 

disturbances in the modal domain; C is partitioned into and C R 

CM figure sensor estimate of C 

CN N x 1 vector which contains the elements of C which are being controlled 

CN figure sensor estimate of the N modal coefficients corresponding to the 

controlled modes 

CN steady -state or final value of C^ 1 

So 

C R R x 1 vector which contains the remaining elements of C 

C* steady-state or final value of C R 

D n N x N diagonal matrix which contains the control system compensation (also 

referred to as the diagonal controller) 

E performance index under the assumption of uncorrelated errors 


( 


Ejf E for a particular set of N actuators 

f frequency 

H M v N matrix which converts the actuator forces to modal coefficients 

retaining the dimensions of force; H is partitioned into H N and H 1 * 

H 1 * N x N matrix which contains the rows of the H matrix corresponding to 

the N modes being controlled 

RxN matrix containing the remaining elements of the H matrix 
hy ijth element of the H matrix 

performance indices 
K gain constant 

M total number of modes (eigenvectors) used to model the mirror 

m diagonal matrix of elemental masses Amj 

m t total mass of mirror 

Am^ mass of ith element of the structural model 

N number of actuators or number of controlled modes 

P total force distribution on the mirror from N actuators 

q vector of disturbance coefficients in the modal domain 

q^ modal error coefficient of the ith mode 

qk* the M x 1 vector q 

qN N x 1 vector of disturbances in the modal domain which correspond to the 

controlled modes 
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R x 1 vector of disturbances which remain after q is partitioned into qN 
and q^ 

remaining modes, R = M - N 

Laplace operator 

M x M matrix of eigenvectors 

M x 1 vector of mirror errors at the grid points 

figure error of the mirror at the ith grid point 

x,y coordinates on mirror at time t 

coordinate axis directed (positive) along the optical axis 

N x 1 vector of forces applied to the mirror surface ( a = arN) 

force distribution over the pad area of the jth actuator 

area over which the pads act 

structural damping of the ith mode 

M x 1 vector of eigenvalues 

wavelength of light 

figure sensor error in determining C M 

density of the mirror expressed in terms of its area 

variance of the ith modal error 

time constants 

defined by equation (59); under the assumption of uncorrelated errors, this 
quantity gives the fraction of '.ie ith modal error which appears in the final 
error 
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^?N for a particular set of N actuators 

& R x N matrix A r H R [a N H N ] 

u> natural frequency 

Subscripts: 


i general term of a vector j 

i,j general element of a matrix 

M last calculated mode ! 

N last mode or actuator under consideration i 

n nth term of a set 

Superscripts: 

- estimate 


T transpose of a matrix 


CONTROL SYSTEM DESCRIPTION 


The Modal Control Concept 

The modal control concept, a applied to mirrors for use in orbiting telescopes, is 
treated in detail in reference 3, and design examples for flat plates are presented. For 
purposes of analysis in the present study, the mirror is considered to be a structure tied 
to a set of supports or mounts that prevent rigid-body motions. The elasticity of the 
mounts themselves may or may not be considered, depending upon the degree of sophisti- 
cation of *he analysis. (The analysis used throughout this paper considers the mounts to 
be rigid.) The modes oi vibration of the mirror, subject to the constraints of the supports, 
are the modes used in the analysis. The mode shapes are referred to as eigenvectors of 
the mirror and the frequencies are the eigenvalues. 


Generally the mode shapes and frequencies must be obtained by numerical methods 
since the solution of the governing partial differential equation is not available. The eigen- 
vectors and eigenvalues oi the structure that are used in this paper were obtained from a 
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numerical program (SAMIS, ref. 4) and have been checked by NASTRAN (ref. 5). These 
results have been verified experimentally (ref. 6). The eigenvectors obtained from the 
S numerical program are tabulated in che U matrix (see appendix B), with column 1 

denoting the first, or lowest frequency, eigenvector and each succeeding column denoting 
. higher order mode shapes. The vector of eigenvalues A (see appendix D) is ordered 
with the lowest frequency first. 

The finite-element model that was used in SAMS is given in reference 7. This 
model was used to extract the first 58 eigenvectors and eigenvalues of the mirror. (See 
appendix D.) This set of eigenvectors and eigenvalues has been used throughout the 
analysis. 

One motivation for using the modal control law was to allow the designer to decouple 
the dynamic behavior of the control system; another and more important aspect of this 
control law is that the mode shapes provide a hierarchy of errors that are likely to occur 
in practice. That is, the modes may be ordered in such a way that mode 1, or the funda- 
mental mode shape, is more likely to occur than mode 2. Also, a measure of the relative 
amplitude is available by examining the eigenvalues of the two modes. That is, if the 
eigenvalues of modes 1 and 2 differ by a factor ox n, the second mode will require about 
n times the input force disturbance to produce the same displacement error. This is 
just another way of saying that the mirror (plant) acts as a filter to high-order modes. 

The one exception to this is that careless initial polishing and figuring of the mirror could 
generate considerable error (as displacements) in the high-order modes. Conversely, 
this knowledge of the mirror should be used to avoid fabrication errors which will be par- 
ticularly difficult to correct. 

System Configuration 

For the purposes of designing a control system for a mirror, the designer obtains a 
transformation from mirror surface deflections (or errors) to modal coefficients, which 
can be viewed as a coordinate transformation. That is, 

q=[u] _1 W (1) 

represents a transformation from the error at a set of points W over the surface of the 
mirror to a ^et of modal coordinates q. Figure 1 shows a block diagram of the mirror, 
figure error sensor, and actuators as they appear in a finite modal representation. The 
mirror itself is mathematically represented by the five blocks (matrices) labeled 
and H^, and A^, and U^. The superscripts N and R have the relationship 

N + R = M 
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where 

M number of modes used to model the mirror (although there are an infinite 

number of mirror modes, practical limitations require a finite number, 
and 58 will be used later for numerical evaluation) 

N number of actuators used (numerically equal to the number of controlled 

modes) 

R remaining modes 

in the physical world the actuator forces are translated directly into nirror 

figure displacements W(x,y,t); in too mnihemai’cal mc<iel the N forces arc transformed 
by the a id H R matrices into a set of fore-, coefficients and a R , respectively, 

in the modal domain. These forces are then transtormed by the A matrices into the 
modal coefficients of displacement. These coefficients, generated by the control system, 
are summed with the coefficients representing the error in the modal domain that previ- 
ously txisted on the mirror, and the result (denoted by and C**) is *r«msformed by 

[UM] into the final displacement W(x,y,t) according to the relationship 



= W M 


( 2 } 


To combine this model into a control system requires a sensor to measure 
W(x.y,t). The sensor output is then changed into the modal coordinate system by the 

proper transformation [u . Since N actuators can control ouiy N modes, the 
subset of the N selected modes to be controlled is usually ail thai is ^e-’erated. One 
would normally control only the first, or lowest, N modes. 


The N selected modes are then fed through the dynamic compensation D^(s) In 
A'hich the proper gains and compensation are applied to each mode independently. If a 
type I system is used, as will be specified in the section entitled "Evaluation of Steady- 
State Errors," then each diegonal element of D 5 ' 1 ’ corresponding to one channel of the 
decoupled controller will contain an integration. The output of D^, denoted a*'*, Is still 
in the* modal domain. In fact, this output is a set of modal coefficients which describe the 
desired force patterns to be distributed on the mirror. To change these to discrete 
forces, which is the way they must be applied to the mirror, the values of a^ must be 


transformed by multiplying by 



This matrix 



also accounts for the effect 



\ 


t 


0 


of the physical mechanism through which the actuator applies a load to the plant. This 
completes the description of the control system which will be analyzed in later sections. 
For a more complete and rigccr>us discussion, see reference 3. 


Modal Representation of Figure Error 

The modal control technique has been developed for mirrors whose modes, denoted 
by Ui(x,y), are assumed to be mercbern of a complete orthonormal set. Since this set 
of modes is complete, a modal expansion of any shape *he mirror can take can be obtained 
by using the correct coefficient for each mode: 

W(x,y,t) * £ CjttWifx.y} (3) 

1*1 


Therefore, for any N points on the minor, li ? fuiiowiitg equation may be written .sym- 
bolically in matrix format: 

r ,i r ir i r .~ 

W^i.yi.t) Ul(*l.yi) • • u s' x l-> f l) CilDj I > C ; (t,Ui Xj.yj':! 

, Cj{t) I 1 


W'x N ,y N .t) U 1 (x N ,y N ) . . . U N > N ,y N ,jjc N W j > C 1 (t)U i (x s .> N , | 

L - L JL j 1 J 


w N - u N c N ♦ u R c R 


where 



• 

u m( x 1 *y j; 

•ir .1 
• # + 

* f 
i 1 

u R c R * 

• 

• 

U Nfl /x N ,3, 'N} • * • 

t'M'xK.yv 

i 

. . i 

• 

i 

j 


and U R is an N x « matrix and C R is an vector. 


{€> 
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The second term in equation (5) causes an error when the mode amplitudes are 
determined since they are not available, and the amplitudes C N are estimates denoted 

CN 

C N = [u N ] 1 W K = C N 1 - [u N ] 1 U R C R (7) 

The term [u N ] U R C R represents the error in the estimate of C. It is possible to 
take more measurements than the number of actuators used. If for example the number 
of measurements is selected to be M (M > N), then 


where 


C't = [uM] 1 W M = C M + 4 m 


4 M = 

u M+i(*i»yi) • • • 
• 

c m+iM 

• 


• 

UM+l( x N>yN) • • • 

• 


( 8 ) 


(9) 


If M is sufficiently large, then will be negligible. Therefore, it will be assumed 

that 


- o 

From figure 1, the following equation may be written (with the x,y,t notation 
dropped): 


C N = q N + A N H N o N (10) 

Since 

a N = -[h n ] D n C n (11) 

and it is assumed that 

C N = C N (12) 
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substituting equations (11) and (12) into equation (10) gives 

C N = q N - A N H N [H N ] X dNc N (13) 

or 

C N = [l -i- A n D n ] q N (14) 

From the other path in the system model in figure 1, 

C R = .\ r H r oN + q R (15) 

Substituting equations (11) and (12) into equation (15) yields 

C R = -A r H r [h n ] f q R (16) 

Substituting equation (14) into equation (16) to eliminate C^ 1 gives 

C R = -A R H r [h N ] + A N D N j V* + q R (17) 

By inserting [A N j \\ N into equation (17), the following expression for C R is obtained: 

C R = -A r K R [h n ]’ 1 [a N ]” 1 A n D N [i + A N D N ] _1 q N + q R (18) 


Equations (14) and (18) give the dynamic values of the modal coefficients of the error in 
the mirror surface. In the present application it Is anticipated that the primary errors 
will be the initial figuring errors and thermal gradients that vary relatively slowly with 
time. Therefore it is reasonable to expect that the system will be generally at or near 
its steady state. The steady-staie performance of the system is discussed in the follow- 
ing section. 

Evaluation of Steady-State Errors 

In determining the steady-state performance of the entire system, first equation (14) 
will be used to assess the resulting error in the controlled modes and then equation (18) 
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will be used to determine the error in the remaining modes. Taking the Laplace trans- 
form of equation (14) and considering the disturbance vector q as a step input allows 
application of final-value theorem to determine the steady-state condition: 


= am c N (t) = lira sC N (s) = lim s[l + A N D N ] colf^j 
ss t-~ s— 0 s— 0 L 1 l s J 


The matrices A* 1 and D** are both diagonal with 


(1 = 1,..., N) (19) 


A n = diag 


O 9 

s' 5 + qs + 


(i = 1, . . ., N; e > 0) (20) 


where equation (20) assumes some structural damping. 

The diagonal matrix D N can be formulated at the discretion of the designer; how- 
ever, a type 1 system is assumed, so that the combination of A N and D N is of the 
form 


A n D n = diag 


K(V + l) . . . (r n s + l) 
s ( r l' s + i) • • • [ r n' s + i) 


n n 

r- S ”| T S + l) + K | t -f 

[l , aNdn] . diag -1 != l± : 

S f [( T ‘' 5 + 1 ) 


which can be simplified to the form 


K |T (vs + l) 

! I + A N D^ = diag — 

s TT [r/s + l) 

i=l ' 1 ‘ 

_ - 
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by properly combining the numerator of equation (22) and factoring. Putting equation (23) 
into equation (19) and taking account of the inverse yields 


CN = ii m /s diag 
ss s-0 ) 


H( Ti ' s + *) 

n-f 1 

n^.i) 


col* 

1 S i 


(i = 1 N) (24) 


This is the expected result that a type 1 system will drive the error in the controlled 
modes to zero. 

In anticipation of evaluating the steady state of equation (18), equation (21) and the 
inverse of equation (23) are combined to get 


r l-l kFT./'tjs + 1) s n (Vs , l) 

lim A N D N |l 4- A N D N = Um diag — 1=1 

s-0 L J s-0 JL. n±i 

silk's + l) kTT (rf s + l) 


Furthermore, since 


lim A® = lim diag — - = diag -i- 

E— 0 s— 0 + £^s + wf 


lim [a n ] = lim diag 
s-0 L -> s-0 


— 0 S^ + CiS + 


= diag 


equation (25) can be used to determine the steady-state value of equation (18): 


C R = lim sC^j = -A r H R |h n ] 1 r.\ N ]’ I q N + q R 
s— 0 J 

C R = q R - A R H R [\ N H N j' 1 q N 
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where equations (26) and (27) indicate the nature of and [a^] . The nature of the 

H matrix and how to evaluate it is given in the next section. Equation (29) states that 
the final steady-state error consists of two parts. The first part (q R ) is that due to the 
original error in the R = M - N modes which were not controlled. The second part of 
the error is that generated by the control system itself as it corrects the error in the 
first N modes. 

A few notes on the dimensionality of the matrices in equation (29) are in order. 

The error vector q N is the initial error in the N modes (not necessarily the first N) 
selected to be controlled and is Nxl; qR is the set of errors in the remaining modes 
and is R x 1. Therefore, 




The matrix is an N x N diagonal matrix which consists of the natural frequencies 

of the modes being controlled; A® is an (M - N) x (M - N) diagonal matrix of the fre- 
quencies of the remaining modes. The and matrices are N x N and 

(M - N) x N, respectively. How these are obtained is given in the next section. 

Determination of Pad Effects 

The function of the H matrices (h n and HR) is to take the point loads of the 
actuators and transform them into modal coordinates. To determine the elements of the 
two H matrices, consider first the continuous case. Let the force distribution on the 
mirror P(x,y,t) be denoted by P 

N 

P = £ aj (t)6j(x,y) (31) 

j=l 

where <*j(t) is the time-varying coefficient of the jth actuator and /3j(x,yl is the dis- 
tribution of the force on the mirror because of the pad. A "pad” is the physical device 
that connects the actuator to the mirror. This force distribution may be expanded by 
using the complete orthonormal set of modes: 


<!, = V hyU, 
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where the (x,y,t) notation is understood. Because of the properties of this set, the coef- 
ficients may be determined immediately: 


»ii-I r Wi dr 


(33) 


where T is the area over which the pad acts. 

Substituting equation (32) into equation (31) gives 


^ v v/’v 
p ■ Z “J 2 h « u ‘ ■ Z Z 

j=l i=l i=l y=l 

Another way of writing equation (31) is 


hijQtj Ui 


(34) 


P = £ a t Ui 
i=l 


(35) 


which expresses P directly in terms of the eigenvectors. The implication of equa- 
tions (35) and (34) is that 


a.- = 


j=l 


h ij a j 


(36) 


or 


fa) = M .foil r N ) 

IJ l a R j H R l i 


\ 


The elements of the H matrix are defined by equation (33). If small pads are assumed 
and it is further assumed that the load distribution is uniform, then 


Pi 


' I, 


(37) 


dA 
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The assumption of small pads leads to the conclusion that the mode shape is relatively 
constant over the area of one pad, and equation (33) becomes 


and 


>H)-Lw r = eij>i<* r - 




ftA 


Ul^l-yi) u l( x 2>>'2) 
u 2( x l*>'l) U 2( X2,y2 ) 


lim H =- 

r-o 


u i( x N’>'n) 

u 2( x N’}’n) 


U M( x l-yi) u M( x 2.y2) • • • u M( x N»yNj 


(38) 


(39) 


When selecting the terms Uj^Xj.yj^ in equation (39). the coordinates (xj.yjj are drier- 
mined by the actuator locations since j3j is assumed zero everywhere except at the actu- 
ator location. Note that the H matrix is nor. .quart. There will be N columns corre- 
sponding to the N actuator locations used; however, all M rows will be present since 
each and every actuator will, in general, excite all M modes. For analysis purposes it 
is convenient to partition the H matrix into two parts; 


H = 


(40) 


The first matrix H R consists of the N rows which correspond to the modes which have 
been selected to be controlled by the N actuators; H R is therefore square. It is not 
necessary that *he first N rovs (N modes) be selected; however, this is usually desired. 
The reason for this will become clear in the examples. This arrangement will be assumed 
in future notation for the H ^ matrix. 

Since H N , H R , A N , A R , q N , and q R have now been obtained, the steady-state 
enor from equation (29) may be calculated for a given choice of actuator locations and 
modes to be controlled. The only remaining consideration is the performance index. 
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PERFORMANCE EVALUATION 


Calculation of the Performance Index 

To obtain the best performance from an optical system, it is necessary to minimize 
the rn.s surface error of the elements (ref. 8). This error is defined as 



For analysis purposes it is usually easier to work with a slightly different quantity which 
provides an equally valid measure of relative performance: 

J = J* 2 (42) 

Using J ai the performance index and changing to discrete notation because of the 
numerical nature of the mirror problem gives 


i=l 

or, in matrix notation, 

Jaiw T [4A]W 

A 


(43) 


(44) 


where 

\ 

[AA] diag [AAiJ 

. J 

W M x 1 vector of mirror displacements or errors 

A total area 

The first N modes will contribute no steady-state error to a step response in a , 

type 1 system since it was assumed that C^ » C^. The final value of W is therefore 
gi*'en by 


18 | • 




(45) 


n 

where Cg S is the (M - N) X 1 or R X 1 vector of the final error in the uncontrolled 
model and is obtained from equation (29). H U, the matrix cf eigenvectors, is obtained 
from a finite-element program — as was done herein — then the eigenvector matrix is 
orthogonal with respect to the mass matrix m (refs. 4 and 5): 


U T mU = I 


For a homogeneous mirror of uniform thickness, m may be specified as an area asso- 
ciated with each grid point in the analysis times an area density constant p. Then m 
may be written as 

m = p diag|AAi] = p[AA] (46) 


or 


i[m].[AA] (47) 

Using equations (47) and (45) in equation (44) giv:s 

J = i^ S ] T U T (48) 

J =^^ uT § U ( C ss] <«> 

j ”^T c «] (6o) 

J - = x wT M w < 5i > 


\ 
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where the total mass is 
M 

m t = ^ Anij (52) 

i=l 

The contribution to the mean-square error of any mode is seen from equation (51) to be 
given by C^m t . 

In any approach where the controller is designed by use of an alternate representa- 
tion of the mirror, it is important to express the desired system performance within the 
framework of the alternate reference frame. The significance of equation (51) is that it 
expresses the figure of merit of the system performance — rms error — as a very simple 
function of the amplitudes of the higher order modes. Thus, minimum rms surface error 
on the mirror is obtained by minimizing the sum of the squares of the amplitudes of the 
higher order modes. If only a relative measure of one actuator arrangement over another 
is of interest, the term l/nit may also be dropped since it is constant for any given mir- 
ror. This gives 



Treatment of Initial Errors or Disturbances 

The q vector obtained from equal, jn (1) assumes that the error at each grid point 
on the mirror surface is known. When this is used in equation (29), the resulting design 
represents a deterministic treatment of the errors. While such a treatment of the errors 
will lead to minimum final error, the result can also lead to ovoroptimism on the part of 
the designer. Consider the following case. 

Given a set of initial errors, the designer determines that a specific actuator \ 

ai -angement will reduce the final error to an acceptable value. When the mirror is 
placed in orbit, it is highly likely that the errors will be different from those anticipated. 

As a result, the second tern in equation (29) — the error generated by the control 
system - will change, possibly significantly. Consequently, the total error as given by 
equation (29) may now be unacceptable. It is concluded, therefore, that the "best 1 ’ actua- 
tor location, based on deterministic errors, is sensitive to initial error. The question is, 
of course, how sensitive? Usually very sensitive, because the actuator positions have 
been chosen to generate specific amounts of error in the uncontrolled modes (generally? 
opposite and equal to what was originally there) when specific amounts of error are 
removed in the controlled modes. A slight change in the error in the controlled modes, 
therefore, could make a great deal of difference in the final error. 
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An alternate way to treat an error is in an uncorrelated fashion, as suggested in 
reference3. In tills approach the values of C^ s are still given by equation (29); how- 
ever, the performance index is now the expected value of the mean-square error 


E(J) 


1_ 

A 


e|w t aa\,vJ 


where E is the expectation operator. If the errors are uncorrelated, then 


(54) 


E(q iqj ) = 0 

and the performance index becomes 


(i * j) 


(55) 


E(J) = E^[q R ] T qR + [^^1^ 
For the assumed type 1 system, where 

* » a r h b (a n h h ]’ 1 


(56) 


(57) 


equation (56) may be rewritten as 


E(J) = 


I 


i=N+l 



M-N 

y i 




l 

i=i 



(58) 


where the variance of the error, which is assumed to have zero mean, is given by a £ . 

Hi 

By reversing the order of summation in the second term of equation (58) and making tire 
additional substitution 


\ 



M-N 

V 



(59) 


equation (58) can be rewritten 
M 


« & 

= X °q, + £j 

i=N+l i=l 


(60) 

Ol 


One notation addition is needed in equation (60), that is, to add the subscript N to 
E and to to indicate the number of actuators being used. This will prevent confu- 

sion later. Equation (60) is then written 


E N = I ^ + Z (6 > 

i=N+l 1 i=l 1 

It should be pointed out that equation (61) represents an expected error. Depending upon 
the inclination of a particular individual, he may choose J*, J, or given by the pre- 
vious equations, as a measure of the performance. The only problem with these equations 
for the performance is that they require exact knowledge of the error vector and may be 
quite sensitive to changes in the error vector, while equation (61) requires only a knowl- 
edge of the variance of the error. It is more realistic to mnke an engineering estimate of 
this latter qiantity than of the actual errors. 

In equation (61) both parts of the error are seen to be positive. This means that the 
two error components will add directly. To influence the expected error, the designer 
may do two things. First, ha should encourage the opticians to keep the error in the 
higher order modes as small as possible because he cannot do anything to reduce this 
error except possibly increase the value of N. Second, he should select actuator loca- 
tions which would minimize the value of In fact, if could be made zero for 

i = 1, . . ., N, then the expected error would be independent of the initial error in the con- 
trolled imdes. Since most of the error will likely occur in the first N modes, choosing 
actuator locations to minimize will lead to locations which tend to produce perfor- 

mance indices virtually independent of changes in error. 

In the following sections and the appendixes, various design examples will be given 
which are based upon the theory developed up to this point. Effects of initial errors, actu- 
ator placement, error treatment, and the number of actuators necessary for a particular 
case will be discussed. 

DESIGN PROCEDURE 


Numerical and Physical Data 

The mirror which will be used in the analysis is shown in figure 2. It is a 1.27-cm- 
thick (’/2-inch), 76-cm-diameter (30-inch), F/3 spherical mirror which is supported on 
a kinematic (non-overconstralned) mount. 

As mentioned in the Introduction, it is generally not possible to obtain closed-form 
expressions of the eigenfunctions of a practical mirror configuration such as that con- 
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altered her?. Therefore, estimates of the eigenfunctions are sought :hrough numerical 
echniquc a. An analysis of the mirror was made with the SAMIS structural ae.ai.vaia pro- 
gram (ref. 4 ) and the grid breakup shown in figure 3(a). This grid breakup wrs adopted to 
comply with the SAMIS recommendation for use of equilateral triangles for maximum 
accuracy. This pattern also insured the location of a grid point (where forces may he 
applied and displacements censed) at each actuator and sensor location of ihe existing 
mirror shown in figure 2. The actuator and sensor locations are shown tn figure 3(b), 
where the numbers refer to grid {joints. As reported in reference 7, tins procedure pro- 
vider, a satisfactory discrete model of th«- mirror behavior. 

In *hc preceding sections, M has represented the finite number r.f modes used m 
lieu of the infinite number which would tie required for an exact representation of the mir- 
ror. In {'"•-forming a specific design it is desirable to select i value for M which is 
high enough to ltirut the 1 error in the .mulvMS and yet is not so high mat it would require 
an exet nsive amount of calculation. Because oi tne inability to place a bound on the error 
incurred through the choice of a specific value of M. the choice is not subject to precise 
evaluation, in the present study me mirror shown m figure 2 was available for corrobora- 
tion of theoretical Jesuits. Tins mirror has 58 displacement sensors, which similarly 
limit the maximum number of modi amplitudes that can tie estimated. For this reason 
M ~ 58 was selected .'or the present study. It is believed that this nun. her is more than 
adequate lor the present application; however, an exhaustive study on this point was net 
performed and it is i^jsstble that a smaller value of M might also Is* satisfactory. 

The matrix of eigenvectors determined irom this model tt> gtvm tn appendix It, and 
appendix C contains ph>t* of node lines {»«■ a selected of modes. Tie' numbers on these 
plots ;>r«- the valet of the eigenvector* at the grid jaunts. Therefore, etch figure repre- 
sents one- eigenvector. The orientation of the figures in appendix C is the same as that in 
figure 3(t) to allow grid point numbers and eigenvalue* to be correlated. Tt.e eigenvalues 
of the mirror are given in appendix f> and are plotted in figure 4. For completeness 
appendix l) aUo contains a listing of the diagonal mans matrix 

For tlr sign and analysts purposes three error rectors were assumed. The first 
error vector was cbtaimd from reference 2 ami is the error ir figuring the mirror. Tht- 
ftecend and third were artutiarlly generated. The second was generated by assuming a 
set of four forces of 4.44# newtons (i p-oucvti located at end jKHr.ts ui, 20, 21, and 40. 

Tnc third example error was a parabolic error. The modal coefficients qj for these 
three « nor example* are listed In utoe I. 

Control System Design Evaluation 

The basic design procedure consists oi setting up a trial design, e valuating C^. 
from equation (29), calculating ,Ij from equation (53), and comparing this value of Jj 
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with previous or desired results. Ii the results are not satisfactory, then a new design is 
tried. Since many trials will be necessary, the only practical approach is to perform the 
design with the aid of a computer. A program to build the various matrices and to evalu- 
ate eg, and the performance index Jj has been written in FORTRAN JV and is given 
in appendix A. The fiow diagram for the program is given in figure 5, where the steps of 
the design procedure are set out in a straightforward manner. 

As shown in figure 5, the program will first read in values of the eigenvector matrix, 
the eigenvalues, and the initial error, or disturbance, vector. The program must then be 
supplied with the number of actuators N to be used and the placement of these actuators. 
Actuator placement is specified by grid numbers. The selection of actuator location ar.d 
number is the major degree of freedom that the control system designer has, and this 
selection more than any other will influence the value of Jy. The program must then be 
supplied with the number of modes to be controlled (the number of controlled modes must 
equal the number of actuators) and these modes identified. Identification of modes is by 
mode number corresponding to the column numbers of the U matrix. With minor excep- 
tions, these should always be the first N modes. From this point the program will sort 
the A, q, and H matrices and carry* out the calculation of c|* and Jj. The program 
will also do one other task; namely, it will calculate the actuator forces ana final error on 
the basis of making a least-squares fit of the errors to the desired rhape. This allows a 
comparison of the modal control law to an optimal control law for the same actuator loca- 
tions. One restriction on the program, and ultimately the designer’s freedom, is that the 
choice of actuator placement must be restricted to grid points of the structural analysis 
model. Tlis restriction has not been found to be serious in the present r uel. which has 
58 grid points. 

One particular word of caution is in order about this, or any other program, which is 
used to calculate the final error. For a given pad shape and size and a given set of modes 
to be control'ed, selecting the actuator placement pattern fixes the H^ 1 matrix. Since the 
inverse of is part of the controller, the actuator locations must be chosen to insure 

that is nonsingular. A singular (or ill-conditioned) matrix indicates that the 

designer has placed actuators in such a manner that the amplitude of at least one controlled 
mode at these actuator locations is (or is nearly) a linear combination of the amplitudes of 
the other controlled modes. If Is singular, the ucuiators cannot independently control 

the give i modes. If H N is 111 conditioned, the modes can be controlled but only at the 
expense of large applied control forces, which generate considerable error in the higher 
order modes. The designer obviously must avuid these cases; however, spotting a poten- 
tially singular matrix mainly on the 1 asis ot an actuator-placement pattern is almost 
Impossible. The computer especially has ^rouble spotting ihis condition, since for rea- 
sons of numerical roundoff and the particular Inversion procedure used, it will obtain 
"inverses" for very ill-conditioned or even singular matrices. The crogram given in 


Read eigenvectors (U^) 
Read eigenvalues (Aj) 
Read disturbances (q^J 


Select number of 
actuators to use and 
modes to be controlled 


Select the N actuator 
locations 


Form following matrices 


H' 1 -' 

L J .N „R 


, n - x««* [a n h n ] q | 


M 

. , I 

isnil 


Evaluate optimal performance 
a . Qvrt] 1 AH]' ‘ [»H T ]q 
E « q. \H j[\a] T Xl1 j * 1 [' H ] T 11 


irur* r '. - r lev iiagras; of cow -ter program in appendix A. This prograra evaluates the 
prrtcrnttr.ee index for a giver, rirrer and a given, or selected, set of actuator 
locations. 
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appendix A treats this problem by calculating and printing out the normalized determinant 
of H 1 ^. If the normalized determinant is very small relative to 1.0, then is said to 

be ill conditioned. (See ref. 9.) Usually an ill-conditioned H N matrix will result in a 
large value of Jj, which automatically excludes that actuator arrangement. This, how- 
ever, is not always the case. (See, e.g., fig. E10 and associated discussion.) 

Results for Deterministic Errors 

From the firs, set of errors a series of design trials were run with various numbers 
of actuators. For the 58-node mirror model, all possible combinations of one, two, three, 
four, and five actuators were surveyed. Beyond five actuators, the number of possible 
combinations becomes too large to make exhaustive searches. Several design "rules of 
thumb" were tried to choose actuator locations; however, none of these provided values of 
Jq that were considered to be near the minimum in light of the results from the exhaus- 
tive searches carried out for fewer actuators. 

The technique that produced the best results was a gradient-type search which used 
the computer interactively. In this technique an initial actuator piaceme-nt is chosen and 
the output of the computer is presented on a CRT. A perspective view of the mirror is 
also generated which shows the deformed state after control. A series of these perspec- 
tive plots for five actuators is given in figure 6. On the basis of the tabulated data and 
the perspective plot, one actuator is mo^ed one grid point and the program reran to see 
whether a gradient can be set up on Jj to improve !he mirror performance. At most, 
six trials are required to exhaust the possible moves for one actuator. The most impor- 
tant single piece of information turned out to be the perspective plot — especially early in 
the search procedure. This plot allowed the grid point with the largest error tc be easily 
spotted, and the actuator nearest this error was moved. The series in figure 6 took 
approximately 2 hours and improved the performance index Jj from an initial value of 
1136.6 to a value of 100, a factor of 10. Further effort beyond this point failed to provide 
improvement. 

The results obtained from these gradient runs can be compared with the known min- 
imum value = 72, obtained fiom an .xhaustive search. The total number of runs that 
were required in the gradient- c :.irch process was 63. If the program were implemented 
so that the computer made all the choices of actuator placement, the gradient search would 
require about 30 seconds of computer time. The use of the computer in the interactive 
mode, however, allowed considerably mere insight to be gained into what factors were 
affecting performance and what factors were not. 

The results of the exhaustive searches for up to five actuators are given in appen- 
dix E. In each case the best 10 locations are shown. The initial performance index in 
all cases was Jj = 1136.6 and two 'inal errors are giver.. The first is the final error 
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obtained with the modal control law and the second is the result obtained when the actuator 
force is selected to minimize the rms error on the mirror. This is referred to herein as 
the optimal control law. A summary of these resuirs if given in figure 7 for the modal 
control law. This figure shows the miirimum error obtained plotted against the number of 
actuators used. The top curve is for the first error example and the lower curve is for 
the parabolic-orror condition. Note that the 'ertical scale is logarithmic. 

The plot in figure 7 suggests an empirical method of estimating the minimum num- 
ber of actuators that a particular disturbance vector might require for a given mirror and 
performance index. That is, an exhaustive search over the model is made for a limited 
number of actuators and the results extrapolated to the desired performance. This would, 
of course, be equivalent to assuming an exponential decay of error with increasing number 
of actuators. The class of errors and plants for which such an assumption would hold is 
not known. Clearly, the error generated by the four loads of 4.448 newtons (1 pound) at 
the discrete grid points is noc in this class since the final error is zero for four or more 
properly located actuators. 

The problem with the deterministic approach is that first, there is presently no way 
to determine readily the best actuator placement, and second, the placement is very sen- 
sitive to initial error. 


Results for Uncorrelated Errors 

When the designer assumes uncorrelated errors, considerably more can be said 
about where actuators should be placed and how many actuators wiM be needed. 

First, equation (51) is minimized instead of equation (29). Since the first term in 
equation (61) is constant for a particular number of actuators, the goal is to minimize the 
second term, which is the error generated in the uncontrolled modes by the control system. 
Instead of writing a pregram to do this, an alternative procedure which enabled the existing 
computer program to be used was adopted. This procedure yielded results that closely 
approximated those which would be obtained from equation (61). The alternate procedure 
consisted of a deterministic minimization of the second term in equation (29). This cor- 
responds to the intent of the uncorrelated case in eliminating the possibility of a control- 

2 

system-generated error canceling an existing error. If the values of are suffi- 

ciently small, there will only be a trivial difference in performance between the two 
procedures. Sir.ce the second term of equation (61) is to be assumed small, the final 
value of the expected error is the first term of eq uation (61). A pint of this portion of the 
equation is giver in figure 8. 

Since the actuator locations are less sensitive to initial error for this criterion, 
exhaustive searches were run for only the first error example. The results of these 
searches are given in appendix F for one, two, three, and four actuators. The values of 
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ripire ?• - Frrur decay as a function of r.uaber of actuators for dcteriainistlc errors 






for tho best run of each number of actuators are given in table II. Since the values 
of ct>^ are relatively small, the expected error falls right on the predicted values of 
figure 8. 

The choice of possible actuator locations car. be considerably reduced when it is 
desired to minimize the generated error. A look at the actuator location for four actua- 
tors shows that each set lies very close to the node lines of mode 5 and several higher 
modes. As pointed out in reference 3, this means that these inodes will not be excited to 
any groat extent. The first mode that is excited to any great extent by this arrangement 
of four actuators is mode il. The ratio of eigenvalues for mode 11 to mode 4 is 


l ll 4.8? xlCT 4 n 
= '■ ■— — a 0. lot) 


2.95 x 10 


This means that the filtering action by the mirror is about 6 limes as much for mede 11 
as for mode 4, the last controlled mode. 

To test the concept of placing actuators at, or near, nodes of higher order modes 
to minimize the generated error, the case for seven actuators was tried. The seven- 
actuator case was chosen because of two factors. First the eigenvalue plot of figure 4 
shows a distinct jump between modes 7 and 8. This implies good filtering action by the 
mirror at this point. Second, the expected-error plot of figure 8 shows a jump between 
six and seven actuators. This implies that a considerable increase In performance can 
be obtained with seven actuators over sLx actuators. In selecting the actual grid loca- 
tions the node lines for modes 8, 9, and 10 were overlaid and a set of grid points near the 
junction of these modes was selected. {See the last diagram of i»g. F5 in appendix F.) 

This reduced set of grid points was searched for those locations which yielded minimum 
generated error. The results of this search are given in table n /c^mand in me last set 
of figures of appendix F (actuator locations). The expected performance index with the 
best arrangement of seven actuators, obtained by using the minimum generateci-orror cri- 
terion, was Jj = 274.9. Had all the values of been exactly zero, the performance 

index would have been = 250.6. The results, therefore, are very close to the predicted 
values of figure 8. 

It should be noied that this result was obtained by searching only a relative few of 
the possible ?ctuator locations, therefore, It cannot be said that the result represents the 
actual minimum. The Important point is that a relatively short computer run was able 
to establish a set of actuator locations having a final error very close to the absolute 
minimum. 
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Comparison of Modal Control and Optimal Control Lavs 

As was stated earlier the computer program which was used to calculate the final 
errors also calculated the error that would Stave been obtained If an optimal control law 
(least -squares fit to best sphere) had been used. In the modal domain the error is given 
by merging equations (10) and (15) to get 


CN 


q N " 


aN H N 

C R 

s* 

3 r 

4 - 

A R H R 


C * q +• AHor 

11 is desired to find the value of a which minimizes 
Ji « C T C 


Therefore. 


(53) 


(04) 


~ » 2{AH] T {A!{{a + 2(AH] T q (65) 

This yields 

or * 

From equation (66) the value of the forces needed to minimize the value of Jj is 
calculated. The vector C is then calculated from equation (63), and Jj is determined 
from equation (64). The computer program outputs the a vector and the value af Jj. 
This can then be compared with the similar value of Ji unuer the modal control law. 

Numerical comparisons are given in appendixes E and F. Generally, given a set of 
actuator locations that were "good," there was little difference in the final result. The 
counterbalancing features of these two control laws are that the modal control law enable* 1 
the dynamic behavior of the mirror to be considered while the optimal control law would 
never make the mirror worse than It was originally. To create more error than the orig- 
inal error would require a rather gross misplacement of actuators on the part of the 
control-system designer; however, it is possible. 

For example, when the errors of example it in table I are treated deterministically, 
the conclusion is that the best actuator location is in the center of the minor. If now the 
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error distribution changes so that the mode I error and the mode 3 error are inter- 
changed, the final error will be much worse than the original. 


The reason for this is that the actuator was placed on a node of mode 1 while 
instructed to control mode 1. This produces considerable force, which generates con- 
siderable error. In the original distribution it was desirable to generate a lot of extrane- 
ous mode 3 error, but in the modified error distribution It is not desirable to generate 
much mode 3 error. The optimal control iaw, faced with the same situation, would do 
almost nothing in the second case, resulting in almost no change of figure error. For 
either ccnirol iaw the actuator placement was bad. 

It should be recognized that any closed-loop control iaw, except one that has a 
decision-making capability as to relative performance, could result in a mirror figure 
that is worse than the original figure if the actuators are placed in a poor location. The 
ability of the modal control law to "order," or to establish a hierarchy of, likely mirror 
deformations through r«»*nv<»ctors and to provide a measure of relative likelihood 
through the ratio of eigenvalues means that the control-system designer will not place 
actuators in a bad location when using this control law. 

Design Examples 

To be specific, consider the case for four actuators and error example 1. Trial 
actuator locations are grid points 2, 30, 4 2, and 45; and the controlled modes are to be 
modes l to 4 Inclusive. The II N matrix Is therefore given by 


-1.092 x IQ’ 2 
N 5.199 x 10“ 3 
11 * -2.273 x ID' 2 
-1.453 x 10' 2 


9.997 x 10-3 
7.237 x IQ -3 
-4.529 x IQ" 2 
-5.481 x 10“ 2 


-1.307 x 10” 3 
3.454 x 10“2 
-3.302 x 10- 2 
-1.235 x 10-2 


3.009 x 10"3 
-1.016 x 10” 2 ! 
-3.347 X 10’ 2 
-1.229 x 10*2 


(67) 


The following information Is then obtained from the computer program for this case: 

(1) The normalized determinant of li N is 0.23, which is considered satisfactory. 

(2/ The force vector a is col(^2.108, 5.9*0, -3.256, 3.78ljnewton3 
(col£o.47i, 1.33, -0.732, 0.850) lb), given in the s me order as the actuator grid 
locations. 

The remaining important results are tabulated below in both their relative values 
(Jl), as determined ' y the program in appendix A., and their absolute values (J*) as given 
by equation (41). 
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Absolute values 
(rms wavelength. 

Rela tive values A, 0.6328 n m? 

0.64 
0.474 
0.433 

0.354 

0.218 
0.217 

The results correspond to the best possible location for four actuators for the determin- 
istic error criterion and error example 1 . 

For the case )ust given, all factors worked tog • o product- an acceptable design; 
however, it is instructive to consider briefly a case teat does not produce good results. 
For example, it was stated earlier that the first (lowest) N modes were usually con- 
trolled. Now consider a case in which this is not true. In error example 1, table 1, it 
can be seen that the four modes which contain most ot the error arc modes 1, 2. 7, and 10. 
Suppose that the above actuator locations (2, 30, 42, and 45) are used and the control sys- 
tem is designed to drive these four mod.'S to zero instead of the first four. The 
matrix would now tie different. For reference, the first column would now be 

col jj-l.002 x IQ’ 2 0,519 x lO - 2 -0.007 * 10 * 2 5.809 x lO-^ 

This, in Hr elf, presen'i no problems; howover, in equation (62) the ratio of the eigenvalue 
of the first mode which would be excited Ajj to the eigenvalue oi the highest controlled 
mode A 4 was 

An 

— — * 0.1C6 < l 

a 4 

which results m an attenuation of the control -system-generated error. For the cast m 
which the four modes containing the most error are controlled, the ratio 


Total original error 1136.6 

Original error in first four modes ........ 618.0 

Original error in remaining 54 modes 518.0 

Error generated by the control system in the 
54 uncontrolled modes 346.0 

Final error by the modal control law 131 1 

Fin^l error by the optimal control law ...... 130.4 


a 3 3.06 x 10 " 2 
A 10 5.67 x 10-* 


(68> 
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reBults in an amplification of the controi-systera-gcneratcd error. One might therefore 
expect that the final error for controlling modes 1, 2, 7 and 10 could be worse than the 
error for controlling mode* l, 2, 3, and 4. The actual calculated final error for r ontrol- 
hng modes 1, 2, 7, and 10 was found to be 

Jl » 1.83 x 10 5 

which is worse than the original error. It might be argued that a different actuator place- 
ment would improve this answer. Although this is probably true to some extent, a new set 
of actuator locations can only change the H*"* matrix but not the utio in equation (68), 
which is the underlying cause of the prob'em. 

CONCLUDING REMARKS 

A design procedure for selecting actuator locations on thin mirrors wiich arc to bc- 
controlied by a modal control law has been worked out for use with typical nmnevicat data. 
Instructions are given for constructing mathematical models of the system. Two ways oi 
treating ui3turbances are discussed. These two techniques, deterministic and untorre- 
lated, are examined from the standpoint of sensitivity to various mirror errors, deter- 
mining the number cf actuator* required, and means of finding the best locations. For 
the deterministic case it was found that the "beat" actuator locutions (those locutions 
which will mlnlmur the error) are very sensitive to th? error distribution: these loca- 
tion* can presently be found only by exhaustive searches of atl possible actuator locations, 
and the number of actuators required for a specific mirror and specific error can only be 
estimated after much computer time is used. For these reasons it is not r ecu a mended 
that error* be treated deterministically because the exact nature of the h.ul bgore-error 
distribution on the mirror surface will cnange with telescope attitude. 

When the error* are treated a* if they are uncorre luted, the locations are much 
less sensitive to variations in error distribution, an estimate of the number of actmMr* 
required to produce a desired reduction in figure error cun easily t.e made, and locatro is 
which will yield result* near these estimates can he found in a reafonabie manner. At 
present, mis technique is much preferred even though it requires more actuators than the 
deterministic method for a specific assumed error. 

For example, the deterministic case for five actuator* and 58 possible locations 
require* 4.6 x 10® trial*. When the error* are treated as uncorrelated variables, how- 
ever, it l* possible to select a subset of the 58 grid points which will yield results close 
to the theoretical limit. This subset is chosen by selecting potential locations near the 
common node line* of the next few higher modes and making an exhaustive search of these 
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locations. Also, for the uncorrelated treatment it is possible to estimate the number of 
actuators a given mirror and mount configuration will require by using an estimate of 
only the variance of the expected errors. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., November 11, 1972. 
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APPENDIX A 

COMPUTER PROGRAM DESCRIPTION. LISTING. AND PRINTOUT 

This appendix contains a listing of the program to calculate figure error of the mir- 
ror on tiie basis of a type 1 modal controller which controls N modes. The program 
contains a great deal of ' ommenl statements to aid the user, but a few additional comments 
are in order. 

First, one must obtain a set of eigenvalues and the eigenvectors for the particular 
mirror to be analyzed. Tnis Is a major undertaking and should Ik* done with a standard 
structural analysis program. Tno size of the eigenvector matrix is a critical item. Since 
actuators can be placed only at grid points of the structural analysis model, a sufficient 
number of these grid points must be used to allow reasonable flexibility in actuator place- 
ment. Too many grid points will result in a matrix that is too large to handle and requires 
excessive storage. The 58 grid points and the 58 x 58 U matrix "fed in this analysis 
have proved to be reasonable. This results in a storage requirement of 110 000g, which 
may be too large for some systems. 

The next requirement is to e'ltain an error vector. This may be obtained from 
experimental da*a by estimating or determining the mirror error at each grid jioint end 
multiplying by or it may be arttficaily generated, as were examples 2 and 3. !r. 

any case tie program assumes that the error vector is already in modal coordinates. An 
option of multiplying this error vector by a constant (to change units) is provided also. 

The next op ton selects ttie output. The short option is recommended for all runs 
except debugging and examining final runs. An example of the complete output is given 
after the program listing. The short option is the last page of output. 

The final option allow s tie dosigneJ to change a few of the actuators and control 
modes without completely rewriting the data cards. By setting NEW = 1, a single actu- 
ator location (or several locations) can be shifted by one number. For example, if a set 
of seven ac.uators is being run and it is djsired to vary the location of one or two actua- 
tors to several selected points on the mirror, it can be done by u<ing this option. Assume 
that the original data cards contain modes 1, 2, 3, 4, 5, C, and 7 as controlled modes and 
actuator locations of 6, 16, 20, 24, 38, 55, and 50, and it is dc-sired to change an actuator 
location from grid point 50 to grid point 49. In this case the MOCHNGE would be zero and 
LOCHNGE would be 1, and the next card would contain the number 49 in proper format. 

The controlled modes would then be- the same, but the actuator locations would be 6, 16, 

20, 24, 38, 55, and 49. Note that it is necessary to place the actuator(s) to be moved at 
the end of the list in the original data card. 
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APPENDIX A - Continued 

The remainder of th° program sorts matrices, calculates the various quantities 
needed to evaluate the performance, and formats the output. One point near the end of 
the program might cause confusion if data on a different mirror are used, that is, the 
conversion of the force vectors to pounds. The value of 1.25 x 10"® assumed that the 
original errors on the mirror surface were in fringes (X = 0.6228 ftm). If not, the error 
can be scaled in the first part of the program. Another potential trouble source is in the 
use of the eigenvalues. SAMIS eigenvalues are inversely proportional to frequency 
squared, and other programs (NASTRAN) output eigenvalues in frequency directly. Tliis 
can be corrected at line 000467 by changing OMEGA SQ(J) = 10000. /LAMBDA(I) to 
OMEGASQ(J) = 1C000. * (frequency squared). 
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PROGRAM ACTUATE CINPuT, OUTPUT. TAPElOl 
REAL LAM0C A ( 58 ) 

INTEGER OMCUU53). ULCCAT(58) 

DIMENSION U158.58), MU0Cl5d)» LUCS' I 58) * OHEGA>uC5i) » UlNVJel58l. 

I riN(50,58), hR( 58,58), (.NI58), UK<5fcl», HOLD I 3d , so) 

OtHENilGN DUMNY(5d,i.) » IPIVUlI SR) , iNGtX(‘j8,2> 

CIMF.NSICN P«CO (58,5d) , PR002(58). CRS5I58) 

DIMENSION HCLOUPI 58. 3d) 

L’lPENSILN AOIvIOE<58) 

OIMENSIUN ARRArI5B, 38? , vECTOR(5d) 

OlfaNSlLN ALPHAI58), ALP i TUR( 50) 

OIPEhSiCN PHII30). LTdEltd) 

DIMENSION UINVFLXI 5d ) 

M « 0 
N56 * 58 


READ IN TmE 0 (EIGENVECTOR) MATRIX. THE LIGEn VECTJK i ARE COLUMNS 
OF Tm£ MATRIX. 

DO 5 J=>1.58 

READ CIO) lu( I * J ) . I* l .58 ) 

5 CONTINUE 


RE AC IN THE CORRESPONDING EIGENVALUES. 

REAL 203. LAMbCA 
PRINT 20‘» * LAK8CA 

REAU IN A OESG.RIPTlVE HEADING FuK HE tR AUK VECTOR. 

RE AC 2CV, ■ LABEL 
PRINT 280. LABEL 

PESO In ERROR VECTOR ANC THE SCALER MULTIPLltR. 

RE AC *03, UINVOE 
PRINT 205. UINVOE 
RE AO 261. FACTOR 

SCALF THE DEFLECTION VECTOR. 


00 «C3 l»i. 58 


0000 76 
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( 

i 


coot 00 
000102 
000104 


000113 

000121 

000i27 

000131 

000136 


0CC140 
0001 00 
000151 


C001 53 
0 JO 166 
000167 


000171 

CU0204 


oiMvcem « uinvoeiii * factor 

403 CONTINUE 

PRINT 262 » FACTUR, UINVOE 
C 

C READ If. IPRINT. IF IPRINT « o P UNTQUT , • KILL 86 ABBREVIATED TO ONE PAGE 

L PER RUN. OTHERWISE PRINTOUT MILL INCLUDE SEVERAL INTERMEDIATE MATRICES. 

C 

REAO 20 6, IPRINT 
C 

C NUN RUN IS ThE NUhBER CF RUNS ”0 BE HAOE. EACH ONc INVOLVES A 

C DISTINCT SET Or MOOES CCN1 ROLLED AND ACTUATOR LOCATIONS. 

C 

RE AC 206 r NL’MRUNS 
C 

C BEGINNING OF OO-LOOP REICH PROCESSES EACH RUN. 

C 

00 350 I RUNS* i .NUM'VINS 

C 

C RE AC IN NEW. IF NEW * 1, A NEW SET OF NODES A.D ACT. LOCATIONS IS READ IN. 
C 

READ 2C6, NEw 
IF INE» .EQ.l) GO TC 415 
C 

C IF NrH .NE. 1, READ IN THE FOLLOWING TWO VALUES WHICH INDICATE THE NUMBER 
C OF POPES AND ACT. LOCATIONS TO BE CHANGED. 

C 

REAP 2X6, FCCHNGE, LCCHNSK 
IF (PCCHNGE.EO. 0> Cl TO 405 
MOCMNGE * MODE + 1 - MOCHNGE 
C 

C REPIACE THE LAST MOCrtNGE UCCE NUMBERS WITH THE FOLLOWING— 

C 

REAO 206. lOMOOCdl. I»MOCHNGE, MODE! 

*05 IF JLCCHNGE .£0. 0 ) GO TO 420 
LOCHNGE « NACT * 1 - LUCMNGE 
C 

t REPLACE THE LAST LUCHNGE ACT. LOCATIONS WITH THE FOLLOWING — 

C 

REAO 206, (ULOCATII) ,1-LCCHNGE.NACT » 

GC TO 420 
C 

C REAC IN AN ENURE NEW SET CF MOOES AND ACT. LOCATIONS. 

C 


i 

I 

i 

i 


! 

i 

! 

< 

I 

f 
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000205 

000222 

OJ02J6 

000237 

030241 


000243 

0QP244 

000246 


0002 fO 
000251 
000252 
COO 2 54 
000256 
000251 
0 ^ 02 63 
1)00265 
000266 
000271 
00021 3 
000276 
000300 
000302 
000303 
COO J 06 
000310 
000311 
000326 
000341 


415 READ 20 o, MOOE, (ONCCCIII, I-l.MQUE I 

REAC 206. (ULOCATIII, 1«1. MODE t 
NACT * MODE 
420 CO 422 3*1 tMUQfc 

hodci.i ■ oacoc(I) 

LUC AT 111 « QtOCATI II 
422 CONTINUE 

if {pcce .eg. ii go rc <? 

C 

C TriE FCUO^ING COOPS SOAI Trtt MUOES AND ACTUATUR LOCATIONS IN ASCENDING 
L ORDER . 

C 

LIMIT * MODE ~ I 

00 0 I«l» 11MIT 

1 8E GIN * I ♦ I 

00 t J J* Id EG IN. MODE 

IF i:»GOC(II .LT. MOOCCJJI I GO TO 6 

IHULC * N.JDCUI 

MOOCUI * MCCv-UI) 

MOLCUJ) * 1 PUL 0 

6 CCNTIFuE 

JO 7 JJ*ICEGIN, NACT 

IF tLOCArm .LT. LJCATIJJ) I GO TO I 

IHGLU * LOCATI II 

LOC A T ill * COCATUJI 

LOC ATI J Jl « ImOLO 

7 Continue 

8 CONTINUE 

9 IF UFfllNT .EC. 0 1 GC TO 15 
PRINT 201 , MUOE, (MUCC ( J I « J * 1 , MODE I 
PRINT 203, ILUCATU), 1*1, NACT I 

15 MMINLSN* 58 - NACT 
C 

C CONSTRUCT THE HN and kk partitions of the eigenvector matrix, 
c 


000343 

II 

* 1 

000344 

NHN 

« 1 

000346 

NHF 

= 1 

000346 

DO 

4C 1*1,58 

0003 50 

IF 

(I .£ J. MOOC (III) 20, 3C 

000354 

20 II’ 

II ♦ 1 

000356 

IF 

(II .UT . MODE! M0CCUI)*0 

000361 

CQ 

25 3-1, NACT 


i 


I 

I 

I 



X 
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000363 



MOLC « LOCATtH 

000365 



HN tNHN t J) * UlM0U>,II 

000313 


25 

CONTINUE 

C0C3T5 



NHN * NHN ♦ 1 

000376 



GO TC 6C 

000377 


30 

CO 35 J*l, NACT 

000601 



MOLC * LUC AT ( J I 

000603 



HR<NFR,J» * U( MULOi 1 1 

000611 


35 

CONl INUE 

000613 



NrtR * AMR ♦! 

000616 


60 

CONTINUE 

000616 



IF IIPRINT .EQ. 01 GO TO 50 

OC0617 



PRINT 210 

000*23 



00 65 l-l, NACT 

000625 



PRIM 21 1 r (HNIl.JJ, J*l, NACT I 

000661 

r 

65 

CONTINUE 


W 

c 

i* 

NOW, PART IT ION Q BY THE ROUES CONTROLLED. 

00C666 

V 

50 

II * 1 

000665 



NR - 1 

000666 



DO 60 1*1,58 

000650 



If (I .fcQ. MGOCUIJJ 52,56 

000656 


52 

CNUIJ « UINVOEUJ 

000657 



cmegaso ( i 1 1 « ioooo-/LAHacAm 

000661 



ii * n ♦ l 

000662 



GO TO 60 

000663 


56 

GfUNRI « UINVOEIU 

000666 



J « MOCt «• NR 

000667 



OMLGASU ( i I « 10000. /LAM0OA< 11 

000671 



NR « NR+ l 

000673 


60 

CONTINUE 

000675 



IF UPRINT .EQ. 0 1 CO TO 61 

000676 

r 


PRINT 213. ICNill, I*x,NACT» 


V 

c 

r 

STORE MATRIX UN IN HCcO 4N1> THEN INVERT >T. 

000511 


61 

00 65 J-l.MODE 

000513 



00 65 I * l - MCT 

000516 



hOLOII.Jl * HNl I,J» 

000523 


65 

CONTINUE 

000530 



CALL FAT INV CmN, NACT, OUMFY , M, OETFRM, IPIVOT, iNDfcX 




l l SCALE 1 


tn 


s 
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0C0540 


IF Cl SCALE .EC. 01 GC Tu 66 

000541 


PRINT 231 

000545 

r 

GO TC 350 


t 

c 

/» 

PRIMS OUT HN INVERSE. 

000546 

V 

66 IF 1 1 PR I NT .EQ. 0 1 GO TC 675 

00054(7 


PRINT 21» 

000553 


00 67 I«1,NACT 

S005a5 


PRINT 2U. C1INC 1,01,2 » 1, N00E1 

000571 

f 

67 CONTINUE 


€ 

NCRI4AU2E THE CETERFINANT CF HN. 

000574. 

V 

675 CIV * 1. 

000576 


00 7C i-i,«00c 

000577 


A01 VIDE III * 0. 

000600 


CO 6S J*1,NACT 

000602 


AOI VIDE ( 1 1 « AC1VI0EII1 ♦ HCLCC 1 . J 1**2 

000610 


60 CONTINUE 

000612 


AOI VICE ( 1 1 * SORT ( AQIV 1CEI 1 1 1 

000617 


OIV » 01/ * AOI VlOZili 

000621 


10 CONTINUE 

000623 

r 

ONCRNAL » OETERM / OIV 


K 

c 

c 

CHECK TFE INVERSION ROUTINE BY RE-MULTIPLYING 

000625 

CO 72 J>l,MGDE 

00062/ 


00 72 I-i.NACT 

OC0630 


PRCOCI , J) - 0.0 

000633 


oa 72 k-i.kode 

000635 


PR0CCI.J1 * PRuOlI.Jl ♦ MCLOU.Kl * HN(K,J1 

00064.7 


72 CONTINUE 

000656 


IF ( IFRINT .EC. 0 1 GO TO 735 

000657 


PRINT 219 

000663 


00 72 1*1, NACT 

000665 


PRINT 211 > CPR0QI1, Jl.J-1, NODE 1 

000701 

#• 

73 COMINLE 


v 

c 

t 

CAECUEATE THE ALPHA VECTOR UNOER The HCOAE CONTROL law 

000704, 

V* 

735 00 740 1*1, NACT 

000706 


ALPHA! I 1 « 0.0 


I 


5 


j 


\ 
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OCC7C7 

0307X1 

C3072X 

000723 

C0073Q 

000731 

000732 

000734 

COO/37 


CCC741 

030742 

000743 

300/45 

000752 

300757 

OC 3/40 

occ7o; 

00076? 

000771 

000773 

000774 

000777 


0310X4 

0ul0l6 

03X0X7 

00)020 

03X022 

OC1025 

031027 

03X041 

031043 

031046 

00X051 

031055 

031057 


DO 74 J*1,M00E 

ALPHMl) « ALPEA(l) ♦ HlIi) 1 * QNEGASGIJ) * QM3I 
74 CONTINLE 

740 ALPSTCMII « ALFHACll 
AT* * « 3. 

00 140 I»1,NACT 

ATAN * ATAH «■ ALPHA) 1 1 4*2 
14U CONTINUE 

ataa • sofiTi*r*w) 

c 

C CALCULATE THE ALPHA VECTOR UNQEA THE OPTIMAL CONTROL LA™. 

C COMPETE t (LAMMlTR * <LAH*Hl l-l * <LAM*H)TR * Q 
C 

DC 150 1*1. MMlNoSN 

CO 150 J«l, NACT 

MOLD * I ♦ MODE 

HOLCIPCLO.JI * HR i I » J I 
150 CGMJNLE 

DQ 156 Jai*58 
CO 155 )*1 i NACT 

MOLCII-JI « HOLOi l * X 1/ CMEGASQUI 
X55 CONTINUE 

156 CCMINLE 

Xf UPAINT .EC. 0 1 CO TO 157 
PRINT 257 

PRINT 202, ( HOLUt I • J 1 . J*l, NACT » 

C 

C INVENT THE MATRIX, AND NORMAL I2E ITS DETERMINANT la THE CALCULATION 
C Cf ; PM*. 

C 

157 01 V * 1. 

Ob ’69 l» l.NACT 

AC1VH EtI) « 0.0 
CO 167 J« l.NACT 

ARRAy(I.J) » 0.0 I 

CQ 165 M 1,58 5 ( 

A ’’RAY (1,1) * ARRAY (1,3) ♦ HCLOU, II*HClD(K, ti 
165 CCNTINjfc 

AOIVICtU) a AOIVIOEUJ ♦ ARRAT ( I , J 1**2 i 

167 CONTINUE 

AOI V * CE ( X I - S0*7 (ACI VU3et II J ! 

01 V * CIV * ACIVluEU) 

169 CCNfJNvE 


! 


i 

l 


y 
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00 


0J10&1 


c 

c 


OJ1072 

CJ107% 

001075 


001107 

oom? 

001115 

OJllit 


031134 
00113b 
001140 
001142 
001143 
001144 
001147 
001151 
001153 
OJ1154 
COUSo 
001161 
001163 
00 1 1 76 
C012Q0 
001206 
001212 
OJ 1213 
001214 
UJ121T 
001221 


001223 

001224 

001225 


CALL RAT I V I ARRAY. NACT, CORKY, H, OETEKM, IMVoT, INOU, H56» 
l ISCAlE 1 

OtTERHN * CeURH / OlV 
If (IPRINT .1(1. 0 t GO TO 175 
PRINT 247, OETERX, ISCALfc, OETERKK 
C 

C PRINT COT THE INVERSE Of l CLAM*h)T * tLAX"rt» 1 
C 

PRINT 258 
CO 171 1*1, NACT 

PRINT 202, ( ARRAY 1 1, J J , J*1 » NACT 1 

171 CONTINUE 
C 

c 

176 00 18C J* l.MOOfc 

VECTOR (41 * CM 31 
180 CCMINUE 

fCLO « 1 ♦ HCOE 
CO 1S5 J* HOLD, 56 
VECTOR I J) » CM J-HCQEl 
185 CONTINUE 

00 1-S3 I -l.NACT 

ALPHAl 1 1 * 0.0 

CO 193 3*1,58 

PROC 1 1 i J( « O.C 
00 190 K « l.NACT 

PRCCII.Jl • PRC Cl 1 • 11 • ARRAY 1 1 . K 1 * hOLCU.M 
190 CuMIMA 

ALP>/(I» * ALPHA J 1 7 A PHJUll *3) * VECTOR! 3 1 

193 CONTINUE 
ALP4SU' * 0.0 

00 194 1*1 ,NACT 

ALPHJ.LH * ALPhSUH ♦ ALPHtU 1 J**2 

194 CONTINUE 

al slh » scrualphsuhi 

c 

C NO* CAlCULATE t TRANSPOSE L 
C 

DO 155 1*1,58 

VtCTORIll - 0.0 
00 155 3*1 .NACT 


I 

I 


X 


APPENDIX A - Continued 


VECTCMII » VECTOR (II « HO.OU.4I * ALPHAI4I 
195 CONTINUE 

00 156 1*1, MODE 

00125* VECTUfttn - VECTOR IIJ - Chill 

0012*6 196 CONTINUE 

00 1C 50 OC 157 l*l,MMlhuSN 

001252 4 - NODE ♦ 1 

001255 VECTCRUI - VECTOR! 4) - (Mil! 

001256 197 CONTINUE 

0012*0 SUN • 0-0 

001261 DO 156 1-1,59 > 

001262 SUN * SDH ♦ VeCTQRU»**2 ! 

001265 198 CONTINUE ■ 

001266 ETE * 5 CRT I SUM) 

C 

C CALCULATE P<00, THE MATRIX PRODUCT— 

C UAMBCAfRI * H(R I * H(NI INVERSE * LAM80AINI INVERSE. THE MULTIPLICATION 
C IS SimifUO BECAUSE THE LAK8DA MATRICES ARE DIAGONAL. 

C 

001271 00 80 4- 1, MODE 

001272 DC 80 l* *. "MINUSN 

001273 PROClI.H - C.O 

001276 CO J5 N* 1 » NODE 

0C13C0 PROC (I , If - PfOUI.iJ ♦ HNIK.jf 

001312 75 CONTINUE 

00131* P«CDU,4» •* PSOOII,4» • CNEGASO! Jf / OMEGAS41MOOE ♦ I) 

001322 80 CONTINUE 

C 

C- CALCULATE AND PRINT INTERMEDIATE DATA ECR PSI MATRIX. 

C 

001327 if IIPfilNT .EC.- 0 ) CC TO 83 

001330 PRIM ?9* 

001333 00 62 4*1, MODE 

001335 PRINT 299 * 4, <PROOU,.l), 1*1, HrtlNUSN » 

001351 82 CONTINUE 

001356 83 DO 65 4*1, MODE 

0013*0 Piil(4> • 0, 

00136* 00 85 l* l, PMINI/SN 

001363 PrtlU) « PMlUl ♦ PRUOII, Jl**2 

001371 8* CONTINUE 

00137* If I SPRINT .EC. 0 > CO TO 9* 

0013U PRINT 29c , tPrtim,l»t, MODE » 

C 
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031+11 

<m*u 

0014X4 
00 : 4 X 6 
001+35 
001 + X 

on+n 


00X445 

00 XV*T 

001452 


00X454 

001433 

001*56 

0014*1 

0J144J 

001465 

00 X 466 

00X46? 

031472 

001474 

001*76 

octsoo 

OJX5C2 
001303 
OJ1504 
001507 
OOX51X 
00151 J 
UJ1514 
001515 
001320 
0314/2 


C MULTIPLY PROO er THE ERROR VECTOR IN TME CONTROLLED KOOcS. 

C 

94 00 35 I« 1 .MKINLSH 

Pft CC2I ! 1 « 0,0 

OC 35 M l.HODE 

P«CD2I X) * PR002U) ♦ PROOl I.Kl «CN<KI 

95 CCRTINUfc 

l«= ! £ PRINT .CO. 0 1 CO TO 99 

PRIM 230. (PF002UI » 1*1 »KM INUSN 1 
C 

C CALCULATE TME ERROR CF CONTROL a 
C 

9t CC 99 1*1. KKINU5N 

CSS5U f * ClRl 1 1 - PR0O2II1 
99 CONTINUE 

P* 

V 

C CALCULATE TME ORIGINAL ERRGR IN THE CONTROLLED NODES, THE ONCON- 
C TRCttfcC NODES, TME ERROR GENERATE!) BY THE CONTROL SYS TEN. ANO 
C TmE ERROR CF CCNTROL. 

C 

ECN « C. 

oc ica I«x »noo£ 

ECK * ECN * CM1I»»2 

100 CCNTIAtE 

SCRTtCN “ SORT (ECM) 

EUM * 0. 

CO 103 1*1, KNINCSN 

Elk « CUM * QR1U**2 
JOS CCATlRUf 

SORT ELK * SORT lEUH) 

00 • ECK ♦ EUN 
SQRIOD * SORTtOO) 

CRC • 0.0 

OC UO i»l, KKINUSN 
CRC • CRC ♦ PR002;il*»2 
UO CGM1M.E 

SORTCftt * SCRT1CRC1 
CP S « 0.0 

00 115 I • 1 .MMIM.SN 

CRS • C«S ♦ CRSSUI *+2 

m conuale 

SCAURS a SCATICRSI 

c 


f 

s 

I 


If 




1 

i 


\ 




s 

— ■ aah» 
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c expam. tme epko* op ashram *£cm ro 5« elenehts hit* ieaos 
C COP PESPUNaif.G TO T*£ CCNTRC..UJ HuOES. TKEN HUtriPir Or TrtC 0 HATR1X* 

c to cnuouire tk f otfucnaw vector. 

c 

OJ1526 !>Q 120 l»l, **l*uS« 

001525 h OkOlt ,11 « C«SS <i I 

001530 120 CON71M/E 

oc;532 n * i 

©01 533 N * i 

001536 DC 133 1-1.5& 

P01536 If tatCCHU .eo. t» *23, U6 

04 nu in tassm • a.o 

001566 It » 1! * l 

0015*5 CO TC 131 

0015*6 128 CRSilll * r>ClGlN*l» 

001551 h » * ♦ 1 

001552 135 CO.lNut 

00 155* 00 1 3 fa t*l.i8 

OC1556 P*U.2tll - C.O 

00155? Od UC J»l * 56 

001501 P*a'2ll» * P3JC2U1 ♦ Utl.JI * CRSSfJl 

C0151Q 136 CUMlhU 

90,57* If UPAtNT ,£C. 0 » CC rr. uj 

00 1 5 ? 5 f*IM 255, t PS002I U , !•! ,58 1 

C 

C CCf.vtai fO»CE v*C'J£S TC PQUHOl PJR PRtJ*T OUT « 

C 

001 606 137 All* * A I A*’ * i.2>t~09 

0 J i 010 AL^H* » ALPMCK * 1.2 5E-07 

C01611 CO 138 I'i, MCT 

OJ 1013 .HPSTOAU) * ACPSTVdl *■ 1.25E-0? 

001615 *1, P wit II * AcPfAtH • 1.251-09 


00161 7 

130 CONHNV.E 


001621 

PAtNT 

701, 

*00fc, 06JCCU1, 1*1, MOOEJ 

001636 

p*:m 

20 a , 

(tOCAT til, 1*1, NACT » 

001651 

PA INI 

232- 

CNONAl 

001 6 > 7 

Pm I NT 

i* 6 , 

(AlPSTuUll, l«l,N*CT 1 

001672 

If tN*Ct . 

u. i i p«;;»t 27o 

001 701 

pm J NT 

755, 

AT AM 

001 TOT 

P It r. 1 

75 >, 

(ALP'iAUI, l«l, NACT » 

OO.T22 

If l NAC T . 

U. 8 1 PAIM 7 70 

OOilJl 

PMlNl 

75 3, 

AlP^iUN 

ou: 737 

PK I N T 

<51, 

Iff 


y 



w* 
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CJ1775 

umo 

002002 

032002 

OG20J2 

C32O02 
002002 
VO 200 2 

V02CC2 

C02C02 

002002 

vJ 2 o 02 

002002 

002032 

002302 

002002 


002002 

032002 

002002 

032002 

002002 


PRIST 2?5- CO,So«TOo, ELM, SORTECM, tun, SyKltJH, <.kv* S^ATChC* 

1 CRi. SvRTCRS 
}5ti CO*' 1 !*Lf 
SlCF 

202 fORRAT 1/ t8tl5.5 I I 

200 fGft.«At t 6113.3 I 

20* fORMT U~l, ♦ 7h£ £ KtN/ALuE S CURRESPCN01 N., TO ThC EIGENVECTOR 
ImATMa , C?fci3.5/1> 

205 FORMAT UH , // U£i5.*l I 

206 FORMAT (201*1 

2 or : ;»>at u*i, #>.u«aea or ccntaoluo notes ■ *. :3.//v * Twer are 
; mCDES *, 1215* (219*, 1215 1 » 

20 * FORMAT I *h , f ACTUATCR LOCATIONS*. 1213, 1/1$*, U1>1 1 

209 FORMAT TSAlCt 

210 FORMAT UH ,// • HN PARTITION OF THE EIGENVECTOR MATRIX */l 

211 »' CftMAt UH , //{ «£!’>. *| 1 

210 FORMAT at- «7/* Tnh *.N PARTITION Of THE Q VECTOR IS *, // 

X ISU5.*1 1 

218 f-ouHH an . // « /in Inverse USTfO hr ROxS.*/ 1 

219 FORMAT UH ,/(• MATKtx TO CHECK MN INVERSION - SHUottl ttt IDENTITY 

1 matrix. usuc hr R0*S.* / 1 

225 FORMAT liH ,// OTA* *ML AN SJJARE knS LKiiOk.*,/ 3Tx, 

1 -URROft*, /// ax, 'ORIGINAL 01 StUAS ANC£*» *<, 2£li.*, // e,K, 

2 •ORIGINAL EkRuR 1N», Ta, 2Slo.*. 7&x. *Tmc COnTROllEO DOES* It 

5 6». *0R! GlNAL t: RU* IN Tt-t*, 3X, 2£l6.*, t 6A, *IINC ONTFOLIE 

*0 KCCfS* /// *X. «£ajluK OthfcftATEC IN »h£ *, 2Ud.*f f 6X, 

5 * 0 ACt NT Rout to newts dr*, / 6X, aCOnTROC SVSTtM*, // 6A. *ERRCR Q 
6F CChTftUL*, **, 2£1 6. A 1 

230 FOPl-Ai ilh , // ton MSI MATRIX ♦ ERROR VECTOR ( ExRCS VECTOR IS 

1 ! n Thl nocui OCXAIM.I ,f * This CCL- ,m N Y : CTuR 01v£S IhL LR* 
2,/ft CEMPAU3 IN EACH MCOt.* H f I«F15.<»1 » 

231 e ii (* a a t am, ut * Tmc ouerminani appears excessive:* large or 

ISMAIL. Tut PON vAS STjPPtU.* I 

232 fORKAT an ,//U2 .*, ♦ IS ThE VALUE Of ThE hCRMALUEO OETERM IN 

UM o» **■«. IF THIS VAtot IS SMALL * / 152, • CCNPMEO TO 1, THE 

inATsa is UL-cc»€mioNic.* i 

2*6 Format ilh . / * FORCE <r£CTOR O-NOcft TMt MUCAL CfifiTllOL LAH l POOH 

U'SJ.*// (£12.*, TEU.Al 1 

2*7 fuRRAT Uti ,// (U.t, ilh • l3«’i,l’,l'JvJH«l03> IS THE VALUE 

IUF T.t Cl' T£An IN ANT uf Thf. MATRIX ( (LAMi-UT * U.A i*rt) 1 USED IN 
2f INOING ALPHA.,/ £13.5, • IS THE HURMlL 1220 OtTEWiiNANT. If THl 
3 S value IS shall CtMPAfltO TO 1, THE MATRIX IS lLL-e<T«01 TIoNEO.ftl 
250 FORMAT Urf • H * Tn£ FuR(.£ VECTOR UNC£P T Ht OPTIMAL CONTROL LAM 



CJ 2 O 02 
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1 ( PCUNOSI .* » U IS12. 4, 7E15.4) ) 

002002 251 FORMAT (1H » // * TP£ FINAL ERROR UNDER THE OPTIMAL CONTROL LAV 

1 IS *, £15. A » 

802002 253 FORMAT ( 1H , / * THE RSS VALUE OF THE FORCE VECTOR UNDER THE OP 

IT I MAC CONTROL LAM IS*, £15. A. * POUNDS.* ) 

002002 255 FORMAT UH , / * THE RSS VALUE OF THE FORCE VECTOR UNOER THE MU 

10AL CCNTROL LAM IS*, £15.4, * POUNDS.* > 

002002 255 FORMAT UH , U * ThIS COLUMN VECTOR GIVES THE UEFLECTI UN OF TH 

IE MIRROR AT EACH GRID POINT. *, / * UNITS ARE THE SAME AS iHE ORI 

7GINAL V VECTOR.*, // (8EI5.4I > 

002002 257 FORMAT (HI ,// 31H LAMBDA * H LISTED BY ROMS. ,1 

002002 258 FORMAT UH ,// 60H THE INVERSE OF l ILAM*H)T * UAM*HI ). LISTE 

10 BY ROmS. ,/ ) 

002002 260 FORMAT tlHl, 8A10 I 

C02002 261 FGMMAT (F10.4) 

002002 262 FORMAT C1H , III ' ERRuR VECTOR, EACH TERM MULTIPLIED BY *, 

1 5-10.4, // { G 5115. 4} ) 

002002 270 FORMAT (lh I 

002002 294 C URMAT UH , Hill V9P THP MATRIX PRODUCT — LAMBDA (R) * HIM ♦ 

1H(M INVERSE * LAHbGAIM JNVcRSE — LISTED BY COLUMNS. / * THIS 

2 MATRIX IS CALLED THE PSI .~*,TRIX. * i 

002002 256 FORMAT UH , II 60H TEE DIAGONAL ELEMENTS UP Th.- PSI TRANSPOSE 

1* PSI MATRIX. . / * THESE AP.t THE PHI SUCAkEO TERMS.* II 

2 ( 6c 15 .4 ) ) 

002002 299 FORMAT UP , II * COLUMN NO. *, 14, U Il>fc.5.4l » 

002002 ENC 


I 

t 

i 


w 

CO 


y 
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01 

THE EIGENVALUES 

CCRRtSPONDlHG 

TO TEE eigenvector matrix u. 







3.509706-02 

3.5CI80E-02 

3.058606-02 

2.96660E-03 

2. 507606-03 

2.506506-03 

2.697906-03 


‘ 


8.9S600E-Q* 

8.979906-06 

5.66980E-C6 

6.89360E-06 

6.890506-06 

3.66690E-06 

3.663306-06 




i.mooE-o* 

1.796606-06 

1*795106-06 

l. 79100E-06 

1. 703 706-06 

1.360506-06 

1.058006-06 




1.0552 0E-06 

9.19350E-05 

9.18190E-C5 

7 ,697706-05 

7.083201-05 

7.0636O6-05 

5.66260E— 05 




5.29U0E-05 

5. 2IE30E-05 

5.1UC906-05 

5. 09760E-05 

6.503506-05 

3.879706-05 . 

3.878106-05 




• J. 588906-05 

2.35360E-05 

3.367 70E-C5 

3.21660E-05 

2.966806-05 

t .79960E-05 

2. 75 ’,606-05 




2.768006-05 

2.6I600E-C5 

2.230906-05 

2.227606-05 

2.068606-05 

2.06760E-Q5 

2.Q12906-C5 




2.Q11T06-05 

1.95690E-05 

1.869006-05 

1.59650E-05 

1.59520E-05 

1.577706-05 

1.576306-05 




1.311306-05 

1.3C8Z0E-05 







> 


the error vector 

IS UlNV * ERROR ExAPPLE HU. 

1 





13 

13 

H 


-1.71866*01 
3.96296*00 
6.1081E-01 
-1.223 76*00 
2.557CE-01 
3.717TE-01 
-•.55I9E-02 
-2.0575E-01 

-i. 2 iooe*oi 
-1.2562E*0l 
-1.3S58E*00 
-3.A520E-01 
-3, 6355E-C2 
5.6867E-01 
2.6367E-01 
5.1006E-02 

7.6458E«CQ 
9. 23926-Ci. 
-2.06606*00 
-5.0369t-0l 
-7.6660E-C1 
-7.697dE-('I 
6.5CS86-C1 

-1.10096*01 
-9.0766E-01 
-1.1568E*90 
-3.73366*00 
-2.55CTE*00 
5. 1375E-01 
6. 1736E-02 

-3.o?61E*00 
-2.1J 76E-01 
9.78226-02 
3.826ot-01 
7.92U9E-01 
6. 72086-01 
1.69366-01 

2.76126*00 
1.06006 *00 
-2.79516*00 
9.60296-01 
6.03666-01 
6.11766-02 
2.1160E-01 

-1.56986*01 
5.38536*00 
-2.79756-01 
6 .2557E-02 
-8.87686-02 
-3.2221E-01 
-1. 93996-01 

-1.60586*00 
1. 60036-01 
5.33076-02 
5.97666—01 
6.36186-02 
-1.36756-01 
2.69916-01 

0 

K 

> 

i 

o 

§ 



ERROR VECTOR, EACH TERN HUlTlPLiEO BT 1.0000 





r. 

3 

C 

a 



“l.7H6C*01 
3 .952 9E *00 
6.108 16-01 
-1.22376*00 
2.5570E-01 
3.7177E-01 
-8.5519E-02 
-2.0575E-01 

-:.2100E»01 
-1.25626*01 
-1.3S58E*0Q 
-3.6520E-01 
-3.V3556-02 
5.6867E-01 
2.63676-01 
5. 1C06E-02 

7.61586*C0 
9.2292E-C1 
-2.06606*00 
-5.0369t-Cl 
—7 .66&0E-Q1 
— 7.69 78E-01 
6.5098t-0l 

-1.10096*01 
-9.0766t-01 
-1.15686*00 
-3. 73366*00 
-2.55076*00 
5.13956-01 
6. 1736E-02 

-3.o761£*00 
-2. U 766-01 
9.78226-02 
3.82666-01 
7.92096-01 
6.02086-01 
1.69366-01 

2.7612E*00 
1.06006*00 
-2.79816*00 
9.63291 -01 
6. 03666-01 
6.11066-02 
2.1160E-01 

-1.5*98E*01 

5.38536*00 

-2.7975E-01 

6.2557E-02 

-8.87686-02 

-3.222't-Ol 

-1.93996-01 

-1.60586*00 

1.60036-01 

5.33076-02 

5.97666-01 

6.36186-02 

-1.3675t-0i 

2.6991E-01 



> ■ .* 

V ' i ; ^ ' 

* 


y 









NUMBER of ccnt rolled mooes 


TNfcy APE MOOES 
ACTUATOR LOCATIONS 

1 2 
2 30 2 

3 A 5 

l AA A5 



Hfi PARTITION OF TK 

EIGENVECTOR 

MATRIX 



-1.032AE-U2 

9.997 IE-03 

1.3983E-C3 

7. 6825E-0A 

3.0095E-02 

5. 19905-03 

7.2879E-03 

3.0AH6E-02 

-1.1727E-02 

-1.616IE-02 

-2, 2785E-02 

A. 52S8E- 02 

-A.2840E-C8 

-2.2867E-02 

-i.3-*75E-02 

-1.A532E-Q2 

5.A817E-02 

-3.A0EAE-02 

-1.A801E-02 

-1.229UE-0 2 

-5.513 5E-02 

1.8150E-C2 

2.2 A51E-C2 

2.7355E-02 

1.725AE-02 

the cn partition ci 

-1.718aE*01 

TFE 0 VECTOR 
1.2 10UF*C1 

IS 

7.4458£*C0 

-l. 10Q9E*01 

-3.o74lfc*00 

MN INVERSE LISTEC 

er rcms. 




-1.1002E*01 

5 .7 1AI E*00 

— 1.5412E*0l 

6. O853E*0G 

-1.17315*01 

i. 06 oee+oi 

6,3121E*Q0 

2.54646*01 

-3.9150E*Q1 

-2.9JU3E+0U 

- 3. IN** 7E *0C 

2. A583E *C1 

-1.92S2E*C1 

2.0553E+UI 

5. 53 7&6*QO 

-A.3tA3fc*01 

2. 72375*01 

-2.A1 27E+C1 

i-2582E*01 

1.2737E*01 

2.697 IE +01 

A. 2A13E-UI 

-1.2536E*C1 

i.393«E*01 

-3. 8769f *00 
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52 MATRIX TO CHECK MS INVERSION - SHOULD BE IDENTITY MATRIX. USTEO 8Y RJ«S-' 
CD 


1.0000E*00 

-*.77*0E-15 

-5.3291E-15 

S.3291E-15 

-2.66*56.-1 S 

-3.552 7E-I5 

I.CCOOE’CQ 

*.**096-15 

-3.3527E-1S 

-8.8siae-i6 

3.5527E-IS 

1.537/E-1A 

l.G0D0E*OO 

-1.0658E-1* 

-2*t6S5E-15 

0. 

3.E915E-15 

2.13161-14 

1.00006*00 

-6.6613E 16 

5.3291E-15 

-1.2018E-1* 

-8.88186-15 

1.2*3*E-1* 

1.00006*03 

LAMBDA • H LISTtE HY ROnS. 
5.92001E-12 E.*S«62E-11 

2.82AA5E-U 

5.6J*36E-12 

*.366*16-11 



2.19S2C-T3 * IO»*l 0» 1001 IS THE VALUE OF THE DETERMINANT OF THE MATRIX I UAM»M|T * (LAH»H1 I USED IM FINDING ALPHA. 

1.81*606-05 IS THl NORMAL I / ED DETERM IHANT • If THIS VALUE Is SMALL uUhPAREl) TO l. THE MATRIX IS ILL-CONDITIONED. 


TuE INVERSE OF I IlAM*H 1T * ILAM*Hl I. USTEO 9Y RONS. 


2.22*55E»15 
-5.2699*E«1* 
-*.22?*OE«l* 
- 1.6 16*2E *15 
1.0179JC *15 


-5.2699*£*1* 
3,f 7S20EU5 
-2.K272£*15 
-1.22750EU5 
-1.337086*15 


-*.22980E*1* 
-2.102I2E*15 
1. *13*76*15 
1.29t&9E*l5 
*.*7C75E*1* 


-l,«l«*2£»15 
-1.227506*15 
1.29t69E*15 
2. 581356*15 
-*,85P70E+l* 


1.0I7V:»t*15 
-1.3T7U8EH5 
*.*70 75£*t* 
-*.8>Q70tU% 
S.P 16J6E*l* 


THE MATRIX PRODUCT — LAMBOA(R) * M(RI * HINl INVERSE ♦ LAIlbOAlNI INVERSE — USTEO 8T COLOMNS. 
THIS MATRIX IS CAlLEC THE PS! MATRIX. 

COLUMN NO. 1 


•r Mr 
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-*.99166-02 -1. 8712 E- 02 

1.02426-02 1.18*06-03 

- 1 . 41 * ie-ui -5.4800E-04 

3.58136-03 4*13926-03 

2.29936-03 4.62O7t-03 

-9.72426-U4 7.15956-05 

8.42946-04 9. 95406-04 


COLUMN NO. 2 

-6.40956-02 -1.407*6-01 

-8.35156-03 -5.81006-03 

-2.20166-04 3. 74816-03 

-1.152*6-03 1. 59106-03 

*.,15536-04 2.23126-03 

-4.403 86-04 -4.596 7 6-04 

1.72326-04 -1.05876-04 


COLUMN SC. 3 

9.53616-03 -1.65326-01 

3.82986-03 1.26666-02 

-1.64426-03 8.17006-04 

4.22496-03 -1.9 n 4 /L-03 

1.541*6-03 4.59296-03 

9.31936-04 1.00436-03 

-9.72416-05 6. 1 7386-04 


COLUMN NU. 4 

-1.67766-01 1.71286+CO 

-4.02496-32 -1.94576-01 

4.26846-02 —3. 1075E— 02 

-6.15206-02 1.19136-02 

-1.74476-02 -3.51026-02 

-1.13786-02 -1.78246-02 

-i. 33546-03 -5.C560E-03 


COLUMN NO. 5 

1. 11126-01 1.9047E-QI 

-l. 90806.-01 -2.01016-02 

4.25246-03 -4.60436-02 

-1.40116-02 -3.74986-02 

-8.15816-03 -3.18C86-03 

1.79416-02 -3.63096-04 

-8. *6676-03 -5.50436-01 


-2.80606-02 3.30646-02 

6.617’.t-C3 -4.92036-03 

-7.2S12c-03 -6.37056-03 

3.34496-03 3.0714.-04 

1. 35374-03 1.67046-03 

2.68496-04 1.98726-0* 

».;JS9E-C4 1.15256-03 


1.64C96 -02 4.9403E-02 

6.61666-03 -1,30676-03 

-5.35326-01 -4.70386-03 

2. 70**6-01 -1.05716-03 

-1.82', 76-0.1 -2.3106F-05 

1.60596-01 -3.04046-05 

3.904*6-04 -4.44646-04 


-1,15*56-07. 1.97H1E-02 

5.71356-04 8.07436-C4 

-2.33896-03 -4.044*t-03 

3.4 5 746- 03 8 . 46 3 7 6- 04 

1 .07426-03 - 9 60616-04 

-2.44276-03 4.57256-04 

-8.58826-04 1.7 504E-03 


1.77356-01 -2.74556-01 

-5.16C36-C2 -3.03426-02 

-2.29936-C2 1.93266-02 

-1.49546-02 -9.26346-03 

-1.15046-02 1.10436-02 

3.06446-3* -7.90816-03 

1.5 1536 -02 -1.94*1*"02 


-9.72586-02 -4.08«2E-02 

-7.68 1 26 - 02 1.7 2 0 5 6- 02 

3.18e4t-C? -6.30026-04 

-1,86106-02 -3.64196-31 

-8,00766-03 -1.23876-03 

4.30956-03 -7. .'9)16-03 

-2.22856-01 -4. 01126-03 




cn 

oo 


THE OlACOHAt EUREMS Of T(-£ RSI I»*t.SfOS£ • RSI MATRIX 
THESE ARE THE RHl SCUAREO TtMS. 


8 . 32 RAE -03 2 . 8120 E -02 3 . 4231 E -02 3 . 6 * 3 I £*00 I.UWE-OI 


f S I MATRIX • ERROR VECTOR I ERROR VECTjCR IS IK THE MOOAL (THAI*. I 
This COtuHR VECTOR CUES The ERBjR CEhERTTEO IN EACH HOOc. 


3 .* 8 * 5£»00 
1.097 15 *00 
A. 3 W 0 E -01 
7 . 150 * 6-01 
i.iim-m 
8 . 8 ** 86-02 
2 .* 250 £-fl 2 


-1.WS2E581 
2 . 3 * 016*00 
*. 77981-01 
- 9* 193 IE - 02 
3 . 25906-01 
2 . 0 * 3 * 6-01 
4 .C 817 E-Q 2 


-i.*o 3 <.E»oo 
7 . 29526-01 
3.08 * 56 -01 
1 . * 1536-01 
». 628 * 6-01 
-*. 0 . 78 £- 0 l 
- 1 . 8 ** 26-01 


2 . 15 * 0£»00 
3 . 77196-01 
- 7 .* 17 * 6-02 
1 . 2918 E -01 
-l. * 5256-01 
1 . 1 * 60£-01 
2 . 2 * 276-01 


-6* I25E E+00 
*. 7 >*JE -01 
- 1 . * 0576-01 
l. 65176-01 
*.«U 6 - 0 l 
1 . 7 * 996-01 
7 . 36 1 76-02 


THIS COLUMN VECTOR CUES Th£ OEELECTICh Uf T iE MIRROR AT EACH SRIO RUlwT. 
UNITS ARE T«C SAME AS THE CR 16 IHAL M VECTOR. 


- 2 . 23376-01 
*.U 02 f-Ol 
-*. 23896-02 
3 .» 1176-01 
R. 7793 C-Q 2 
- 5 . 75 T 5 C-Q 1 
-l.l**& 6*0 0 
- 1 . 41 * 66-01 


-*.tj: 5 E- 0 l 
- 5 . 8820 C —01 
- 1 . 39226-01 
-*. 0 »* 7 E -01 
- 2 . 73996-02 
-R.SI 73 E -01 
-*.* 85 »E-a 
3 .SCQ 5 E -81 


- 2 . 2 RSRE-C; 
7 . 7151 E-C 2 
2.3703E~vl 
6 . 3 . 1 LE-C 2 
A./OCOfc-Jl 
-2 . 39 1 26 -Cl 
- 2 .SU 2 E -01 


- 2 . 679 * 6-02 
7 . 7709 E -01 
2 . 7 *u 46-02 
2 . 3595 E -01 
-S.ROm-Ol 
-*. 9 i) 20 t-Ql 
- 6 . 09096-02 


-I.WM-OI 

r. S3 m-oi 

1 . 05146*00 
1 . 7534 E -01 
3 .R 3 /SE -01 
2 . * 5 * 46-01 
3 . 59506-01 





1 


1 . 02936 . JO 

- 3 . 17536*00 

- 3 , 926 * 6*00 


i.uase-oi 

1. 90026 *-00 

- 1 . 01216 *00 


8 . 381 / 6-03 

2 . 555 lfc-Ol 

- 2 . 8 * 5 * 6-01 


9 . * 2 151-02 

- 1 . 78 S 5 E -01 
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APPENDIX B 

EIGENVECTORS OF THE MIRROR 

This appendix contains a listing of the U matrix (matrix of eigenvectors) for the 
76-cm-diameter (30-inch) thin mirror. This matrix was obtained from the SAMIS pro- 
gram in reference 4. A more graphic display of the first 10 eigenvectors is contained in 
appendix C. The eigenvalues associated with each eigenvector are given in appendix D. 
Lambda(l) is associated with column (1), and so forth. The diagonal mass matrix m is 
also given in appendix D. The U matrix is orthogonal with respect to the mass matrix 
(refs. 4 and 5): 
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APPENDIX C 

EIGENVECTOR DIAGRAMS 

This appendix contains schematic representations o' the first 50 modes of the thin 
mirror (figs. Ci to CIO;. For each of these modes toe mirror deflection along the Z-axis 
at each grid point Is written beside that grid point. The contour lines indicate the nodes 
of each mode. Note that mod? 3 is a bowl-shaped mode and therefore has no node Une. 
The numbers at the grid points are taken direct iy from appendix B, column* 1 to 10 inclu- 
sive, and are displayed in tide manner to pro.’tde r..*re tns.ght n:o actual shapes than can 
be obtained from looking at columns oi numbers. Each number has been multiplied by iOQ 
for scaling purposes. 














Hg-ure C9.- *Vde 9 of the thin mirrcr 



Figure CIO.- Mode 10 of the thin mirror 
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APPENDIX D 

EIGENVALUES AND MASS MATRIX FOR THE MIRROR 


This appendix contains a listing of the eigenvalues and the diagonal mass matrix m. 
The SAMIS program eigenvalues are inversely proportional to frequency squared and may 
be converted to frequency (in hertz) by the following relationship: 
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The elements of the mass matrix have the units of Ib-secVin. 


The eigenvalues of the U matrix are 
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APPENDIX D - Concluded 


The elements of the diagonal mass matrix are 
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APPENDIX E 

BEST ACTUATOR LOCATIONS - DETERMINISTIC CASE 

This appendix contains an enumeration of actuator locations which were found to 
produce minimum error when the errors were considered to be completely known (deter- 
ministic) and nonvarying (figs. El to E12). Each figure contains 10 diagrams which show 
the best 10 actuator locations for one of the three error examples and a specific number 
of actuators. Grid numbers for the actuator locations can be found by comparison with 
the numbered pattern in the lower right-hand portion of each figure. The values of the 
final error for the modal control law and the optimal control law are given beside each 
figure in the form 


where A is the error under the modal control law and B is the error under the optimal 
control l?w. These final errors are those given by the square root of equation (53), which 
require'; that these values be multiplied by 0.40 to obtain rms error in microinches or by 
0.019 to obtain rms error in wavelengths. 

A particular point of Interest occurred in figure E10 (two actuators, error exam- 
ple 3). In this figure the H N matrix was decidedly ill conditioned for most of the exam- 
ples. Tne normalized determinant was as low as 6 x 10"^ and the best value was 
3.5 x 10*2, jf this case arose in practice, it would be best to look at three actuators or 
more. In figure Ell (three actuators) the normalized determinant was of the order of 0.3, 
which is very good. 
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APPENDIX F 


| 

j 

BEST ACTUATOR LOCATIONS - UNCORRELATED ERRORS 

t This appendix contains an enumeration of actuator locations which were found to 

t generate minimum error under the modal control law (figs. FI to F5). These were 
{ obtained by using the errors of example 1. The answers are given beside each diagram 
is the figures in the following form: 

| A 

’ B 

C 


where 

A the error predicted by equation (61) assuming all values of < are 0 

B the error obtained from equations (53) and (29) 

C the error under the optimal control law 

The performance index is diat obtained from the square root of equation (53) and must be 
multiplied by 0.49 to obtain rms error in mlcroinchec or by 0.019 to obtain rms error in 
wavelengths. 

Selecting other error examples would, of course, result in the selection of different 

actuator locations; however, It can be seen from the values of i n table U that the 

effect of a different actuator location could not make the final answer much better oecause 
2 

the valeee of are already very small. For this reason the searches for actuator 

locations were restricted to the one example. The searches for one to four actuators 
Inclusive considered uii possible combinations, whereas those for seven actuators con- 
sidered only a small subset of all possible combinations. This subset was chosen from 
those locations near the nodes of the next three higher order modes. This reduced the 
number of runs required to a reasonable value and resulted in a selection of actuator loca- 
tions which were reasontoly close to the theoretical limit. 
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,i‘or IceV. i;c-r vhic:. rir.'. r’..-r ► be error prneritrd by the cent -o’ syst < n Sr driving the flr*t 
to In.* one actuator, tnvr distribution is t.aXen from error e>nrplc 1, 
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cure r*.. - Art -a' cr locations vfc!ch result**! in the *inlttws flen*r*tcl error for error ex<t*?1* X. Tb 
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TABLE H.- VALUES OF ^ FOR N - 1. 2, 3, 4, AND 7 

[ ^Determined for error example I; actuator grid locations 

: corresponding to these values are tabulated below the 

values of 
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Actuator grid locations 

52 

34, 21 

16, 23, 47 

16, 20, 24, 55 
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